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The halo of the Miky Way might contain numerous and dense substructures 
inside which the putative weakly interacting massive particles (suggested as the 
main constituent of the astronomical dark matter) would produce a stronger 
, annihilation signal than in the smooth regions. The closer the nearest clump, 

O ■ the larger the positron and antiproton cosmic ray fluxes at the Earth. But the 

I ' actual distribution of these substructures is not known. The predictions on 

O ' the antimatter yields at the Earth are therefore affected by a kind of cosmic 

variance whose analysis is the subject of this contribution. The statistical tools 
■ to achieve that goal are presented and Monte Carlo simulations are compared 

to analytic results. 
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1. Motivations 

The universe contains a significant fraction of invisible and non-baryonic 
matter 1 whose nature is still unresolved. This dark matter could be made of 
weakly interacting massive particles, 2 such as the supersymmetric neutrali- 
nos, whose annihilations inside the galactic halo might produce a sizeable 
yield of positrons and antiprotons. These rare antimatter cosmic rays will be 
detected with improved accuracy by the forthcoming experiments. 3 ' 4 Nu- 
merical simulations of structure formation indicate that dark matter (DM) 
is clumpy on small scales. 5 Because neutralinos are more tightly packed 
inside the DM substructures, their mutual annihilation is enhanced. The 
resulting increase of the associated antimatter signals at the Earth has been 
so far accounted for by an overall boost factor and by shifting upwards the 




Fig. 1. The positron fraction is derived for a 50 GeV Kaluza-Klein inspired particle 6 and compared to the HEAT excess 7 which cannot 
be explained by the pure conventional background 8 in black. An isothermal profile has been assumed for the dark matter halo of the 
Milky Way, a fraction / = 0.2 of which lies in clumps with mass 10 7 Mq and intrinsic boost B c = 200. Two random realizations of that 
clumpy DM halo are featured in blue. The distance of the closest substructure has been set equal to 1 kpc in the left panel and decreased 
to 0.1 kpc in the right panel. The green curve corresponds to the traditional and wrong shift by a factor of / X B c = 40 of the positron 
spectrum which a completely smooth DM halo would yield. 
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cosmic ray fluxes obtained in the case of a smooth distribution of DM par- 
ticles. 

This approach is too simple. If a substructure were to lie in our im- 
mediate vicinity, we would observe a significantly distorted spectrum as 
illustrated in Fig. I for a 50 GeV LZP species. The question then naturally 
arises to know if such a possibility is probable or very rare. Should we know 
the exact location of each DM clump, the positron and antiproton spectra 
would be derived exactly. On the contrary, they are affected by a kind of 
cosmic variance because an infinite number of different halo realizations are 
possible. The theoretical predictions on the detectable fluxes at the Earth 
must take into account that lack of knowledge. 

We have therefore built specific tools 9 to address this issue. We will con- 
centrate here on the pedagogical example where a positron line is produced 
through the mutual annihilations of the DM particles \ m t° e+e~ pairs. 
This scenario is plausible in the framework of Kaluza-Klein inspired theo- 
ries. In the absence of any substructure, the DM distribution p s is smooth 
and the resulting flux is given by 



The probability for a positron injected at x to propagate towards the Earth 
which it reaches with the degraded energy E < Es is described by the 
Green function G(x). In the presence of clumps, the DM distribution is 
given by the superposition p = p' s + Sp where a smooth component p' s is 
still present while Sp accounts for the substructures. The propagator G has 
already been discussed in the literature. 10 Positrons loose energy through 
synchrotron radiation and inverse Compton scattering on stellar light and 
on the CMB. The farther they originate, the smaller their energy E at the 
Earth for a fixed value of Es- Positrons that are detected at the energy 
E have been produced within a typical range Ad which decreases as E 
increases towards Es but is nevertheless much larger than the size of the 
DM substructures. The positron flux becomes 




(1) 



where the factor S is defined as 
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where the contribution from the ith clump is <pi — S x G(xj) x £j. That 
minihalo produces as many positrons as if the entire volume 

{my* (4) 

PO J ith clump I PO J 

were filled with the density p - The volume £j can also be expressed as a 
function of the substructure mass Mj and intrinsic boost Bi. 

An infinite set of halo realizations must be taken into account, each 
of which produces a different flux <f> r . The boost factor B = <j>/(f> a at the 
Earth is not unique and must be treated as a random variable. We present 
in the next section a procedure to determine the statistical law according 
to which B is distributed and will show that its average value and variance 
depend on the energy E. 



2. Computing the Odds of the Galactic Lottery 

The statistical properties of the random flux <p r and of the associated boost 
factor B have been thoroughly investigated in a recent analysis. 9 I will just 
point out its salient features and summarize the hypotheses on which it is 
based. 

(i) To simplify the discussion and without loss of generality, we assume 
identical clumps with mass M c and intrinsic boost B c . The random flux cf> r 
simplifies into 

<j> r = S x £ c x G * > ( 5 ) 

i 

where the effective volume £ c = B c M c /po is the same for each minihalo. 

(ii) The actual distribution of DM substructures is one particular re- 
alization to be taken from the statistical ensemble made up by all the 
possible random distributions. The flux contributed by the clumps and 
the boost factor must be averaged on that infinite set to yield (<fi r ) and 
B c g = (B = (</>/<f> s )) whereas the variance is defined as 

of = ($.) - {4> r ) 2 and <r B = <r r /(f) s . (6) 

(iii) Clumps are distributed independently of each other. This is a strong 
assumption that generally holds because the substructure correlation length 
is smaller than the propagation range Ad- We just need then to determine 
how a single clump is distributed inside the Milky Way halo in order to 
derive the statistical properties of an entire population of Njj such sub- 
structures. If tp denotes the contribution of a single minihalo, we get 

(<f> r ) = N H (p) and of - N H a 1 = N H { {ip 2 ) - (<^) 2 } . (7) 
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Fig. 2. Positrons are injected with the energy Es = 100 GeV and detected at the energy 
E. The effective boost factor B c ff (black line) accounts for the average enhancement of 
the positron signal resulting from DM dumpiness. Although it has been so far considered 
to behave as a constant, the boost actually depends on the energy. Furthermore, as E 
approaches E$ , the boost variance cfb increases significantly. The region from which the 
positrons detected at the Earth originate shrinks and the number N§ of clumps which 
it contains decreases. 



(iv) The set of the random distributions of one single clump inside the 
Milky Way halo T>h makes up the statistical ensemble T which we even- 
tually need to consider. An event from that ensemble consists in a clump 
located at position x within the elementary volume d 3 x and is weighted by 
the probability p(x) d 3 x. The spatial distribution p(x) and the flux density 



February 5, 2008 15:34 



WSPC - Proceedings Trim Size: 9in x 6in salati'rhodes'06 



6 



of probability V (ip) are related by 



V (p>) dip = dP = 



f p(x)d 3 x 



(8) 



where 2? v is the region of the halo in which a clump must lie in order to 
contribute a flux ip at the Earth. We readily conclude that any function T 
of the flux p (for instance p itself or its square <p 2 ) is given on average by 



Once p(x) is known, B e g and ub can be derived. 
3. Results and Perspectives 

We have applied that procedure in the case where p' s = (1 — /) x p s and 
have finally assumed that the clumps trace the smooth DM distribution 
(which is not the case in general) by imposing that p(x) oc p s (x) . The main 
results are summarized in Fig. 2 and 3. 

To commence, the average boost B e s depends on the energy. As E de- 
creases, the part of the halo which contributes effectively to the positron 
signal at the Earth grows. Regions located close to the galactic center come 
into play, with a much larger density p s than in the solar neighborhood. 
At fixed intrinsic boost B c , neutralino annihilation inside clumps is rela- 
tively less enhanced in these dense regions than in our vicinity, hence the 
behaviour exhibited in Fig. 2. 

We also find that the boost variance erg increases significantly as E 
approaches E$- The above mentioned domain of the halo, inside which the 
positrons detected at the Earth with the energy E have been produced, has 
a typical volume Vs — (V2n Ad) 3 which shrinks as E gets near to Eg. The 
number N$ of the substructures which it contains falls down and the boost 
variance increases because 



Finally, in the large N$ regime, the boost distribution follows a 
Maxwellian law as is featured in Fig. 3. This result may be derived with the 
help of the central limit theorem. In the opposite regime where Ng < 1, the 
boost distribution is driven by the product Nh x V (<p) where the density of 
probability V (p) for the flux generated by a single clump comes into play. 

The method outlined here can bc applied in particular to the interme- 
diate mass black hole scenario. 11 The spatial distribution p(x) no longer 




(9) 



<r r ob 



I 



(10) 
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traces the smooth density p s . An effective boost -B e ff of 
found 12 for both antiprotons and positrons with a large variance. 



E e+ = 80 GeV (E s =1 00 GeV) 



140 
120 
100 
80 
60 
40 
20 




boost(80 GeV) 


Entries 


1000 


Mean 


f\ QQQQ 

u.yyyy 


RMS 


0.1101 


X J /ndf 


6.484/16 


K 


39.7S± 1.26 


V 


0.9998 ± O.O035 


| 

o 


0.1096 ±0.0026 



2a n - 



0.2 0.4 0.6 0.8 1 1.2 1.4 

B MC f B eff 



1.6 1.8 2 

(B eff = 20.5) 



Fig. 3. In this Monte— Carlo simulation, clumps contribute a fraction / = 0.2 to the 
mass of the Milky Way DM halo (assumed to follow a NFW profile with a scale radius of 
25 kpc). Each substructure has a mass of 10 5 Mq. A thousand different realizations of the 
DM halo have been generated here, each of them involving 271,488 clumps. The positron 
injection energy is Es = 100 GeV. The number of realizations has been plotted as a 
function of the boost ratio rj = B/B^g. The rms value of 0.1101 yielded by the Monte- 
Carlo for r) is in excellent agreement with the anaytic value a n = cr B /B c f[ = 0.1097. 
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